* use "data_for_Pischedda_Cheon_CPS.dta", clear

 
**********Pischedda and Cheon, "Does plausible deniability work? Assessing the effectiveness of unclaimed coercive acts in the Ukraine war," CSP**********

******Regression tables
***Table 1 (main text)
reg weapons_view unclaimed
estimates store m1, title(Model 1)
reg weapons_view unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m2, title(Model 2)
reg strike_view unclaimed
estimates store m3, title(Model 3) 
reg strike_view unclaimed  ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m4, title(Model 4)
reg war_view unclaimed
estimates store m5, title(Model 5)
reg war_view unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m6, title(Model 6)

esttab m1 m2 m3 m4 m5 m6 using Table1_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted

***Table A1 (online appendix)
reg weapons_view_dummy unclaimed
estimates store m7, title(Model 1)
reg weapons_view_dummy unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m8, title(Model 2)
reg weapons_pref unclaimed
estimates store m9, title(Model 3)
reg weapons_pref unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m10, title(Model 4)
reg weapons_pref_dummy unclaimed
estimates store m11, title(Model 5)
reg weapons_pref_dummy unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m12, title(Model 6)

esttab m7 m8 m9 m10 m11 m12 using TableA1_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted

***Table A2 (online appendix)
reg strike_view_dummy unclaimed
estimates store m13, title(Model 1)
reg strike_view_dummy unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m14, title(Model 2)
reg war_view_dummy unclaimed
estimates store m15, title(Model 3)
reg war_view_dummy unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m16, title(Model 4)

esttab m13 m14 m15 m16 using TableA2_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted

***Table A3 (online appendix)
reg strike_pref unclaimed
estimates store m17, title(Model 1)
reg strike_pref unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m18, title(Model 2)
reg war_pref unclaimed
estimates store m19, title(Model 3)
reg war_pref unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m20, title(Model 4)
reg strike_pref_dummy unclaimed
estimates store m21, title(Model 5)
reg strike_pref_dummy unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m22, title(Model 6)
reg war_pref_dummy unclaimed
estimates store m23, title(Model 7)
reg war_pref_dummy unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin
estimates store m24, title(Model 8)

esttab m17 m18 m19 m20 m21 m22 m23 m24 using TableA3_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted

***Table A4 (online appendix)
reg anger_level unclaimed
estimates store m1, title(Model 1)
reg reputation_scale unclaimed
estimates store m2, title(Model 2)
reg intention_scale unclaimed
estimates store m3, title(Model 3)

esttab m1 m2 m3 using TableA4_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted

***Table A5 (online appendix)
sureg (weapons_view unclaimed ideology_scale int_trust superior force mil_strength going_war trump putin) (strike_view unclaimed ideology_scale shame_level superior force mil_strength going_war trump putin) (war_view unclaimed ideology_scale int_trust shame_level force mil_strength going_war trump putin), corr
estimates store m25, title(Model 25)

esttab m25 using TableA5_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted

***Table A6 (online appendix)
reg weapons_view lives_importance
estimates store m1, title(Model 1)
reg strike_view lives_importance
estimates store m2, title(Model 2)
reg war_view lives_importance
estimates store m3, title(Model 3)

esttab m1 m2 m3 using TableA6_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted

***Table A7 (online appendix)
reg weapons_view unclaimed US
estimates store m1, title(Model 1)
reg weapons_view unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin US
estimates store m2, title(Model 2)
reg strike_view unclaimed US
estimates store m3, title(Model 3)
reg strike_view unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin US
estimates store m4, title(Model 4)
reg war_view unclaimed US
estimates store m5, title(Model 3)
reg war_view unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war trump putin US
estimates store m6, title(Model 4)

esttab m1 m2 m3 m4 m5 m6 using TableA7_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted

***Table A8 (online appendix)
reg weapons_view unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war partisan2 trump putin
estimates store m1, title(Model 1)
reg weapons_view unclaimed ideology_scale int_trust shame_level superior force mil_strength going_war partisan2 
estimates store m2, title(Model 2)

esttab m1 m2 using TableA8_regression.rtf,  cells(b(star fmt(3)) t(par fmt(2))) starlevels(* 0.10 ** 0.05 *** 0.01) noomitted



******Mediation analysis
***Figure 1
medeff (regress attribution unclaimed) (regress weapons_view unclaimed attribution ideology_scale int_trust shame_level superior force mil_strength going_war trump putin), treat(unclaimed) mediate(attribution) sims (1000)

***Figure A1 (online appendix)
medeff (regress attribution unclaimed) (regress weapons_view unclaimed attribution ideology_scale int_trust shame_level superior force mil_strength going_war trump putin death_scale cognition_scale partisan2 trust_gov), treat(unclaimed) mediate(attribution) sims (1000)

***Figure A2 (online appendix)
medeff (regress anger_level unclaimed) (regress weapons_view unclaimed anger_level ideology_scale int_trust shame_level superior force mil_strength going_war trump putin), treat(unclaimed) mediate(anger_level) sims (1000)

***Figure A3 (online appendix)
medeff (regress reputation_scale unclaimed) (regress weapons_view unclaimed reputation_scale ideology_scale int_trust shame_level superior force mil_strength going_war trump putin), treat(unclaimed) mediate(reputation_scale) sims (1000)

***Figure A4 (online appendix)
medeff (regress intention_scale unclaimed) (regress weapons_view unclaimed intention_scale ideology_scale int_trust shame_level superior force mil_strength going_war trump putin), treat(unclaimed) mediate(intention_scale) sims (1000)


